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Abstract — The average weight distribution of a regular low- 
density parity-check (LDPC) code ensemble over a finite field 
is thoroughly analyzed. In particular, a precise asymptotic ap- 
proximation of the average weight distribution is derived for 
the small-weight case, and a series of fundamental qualitative 
properties of the asymptotic growth rate of the average weight 
distribution are proved. Based on this analysis, a general result, 
including all previous results as special cases, is established for 
the minimum distance of individual codes in a regular LDPC 
code ensemble. 

Index Terms — Low-density parity-check (LDPC) codes, mini- 
mum distance, weight distribution. 



I. Introduction 

LOW-DENSITY parity-check (LDPC) codes, originally 
introduced by Gallager [1], are a family of linear codes 
characterized by a sparse parity-check matrix. Owing to 
their capacity-approaching performance under low-complexity 
iterative decoding algorithms, LDPC codes have attracted 
tremendous attention in the past years. To evaluate the the- 
oretical performance of an LDPC code, a typical method is to 
estimate its performance under maximum-likelihood (ML) or 
iterative decoding assumptions. The performance of a linear 
code under ML decoding can be well estimated based on its 
weight distribution 1 1 1, so having the knowledge about weight 
distributions of LDPC codes facilitate the analysis of the ML 
decoding performance. 

The first analysis work on the weight distributions of LDPC 
codes was given by Gallager in his pioneering work [1], where 
he studied the weight distributions of binary regular LDPC 
codes. Moreover, he also generalized the analysis to non- 
binary regular LDPC codes over Z m (m > 2), characterized 
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by zero-one parity-check matrices. Ever since the publication 
of £0, there has been a lot of work extending the analysis 
of weight distributions of binary LDPC codes in different 
ways, such as |2_-|__. A generalization of weight distributions, 
also known as spectra, of regular LDPC codes over finite 
fields and arbitrary abelian groups were later studied in __], 
|__. More recently, the binary weight distributions of non- 
binary LDPC codes also received some attention 1 10]. By now 
a bundle of formulas about weight distributions of various 
LDPC codes is known, but the value and significance of 
most formulas is far from being fully understood, except in 
the case of binary regular LDPC codes, which have been 
well studied _Q, |__. The difficulty is due to the complex 
expressions for the weight distributions of LDPC codes, which 
are usually obtained by the generating function approach and 
hence are typically expressed as coefficients of a polynomial. 
Given a polynomial p(x) with nonnegative coefficients, a usual 
approach for estimating the coefficient of a monomial x k in 
[p{x)] n is to calculate the infimum of [p(x)] n /x k over all 
positive x, which gives an upper bound of the coefficient and in 
fact has the same asymptotic growth rate as the coefficient [__ 
Theorem 1]. However, analyzing functions like inf y> o f(x,y) 
is not an easy job. When f(x,y) is complicated, determining 
the shape, such as monotonicity, convexity, and zeros, of 
inf y> o f(x,y) becomes a difficult mission. 

In this paper, we shall perform such a mission for ensembles 
of regular LDPC codes over finite fields. At first, as an easy 
consequence of the results in |__, __], IfTD . an exact expression 
is introduced for the average weight distribution of a (c, d)- 
regular LDPC code ensemble over the finite field ¥ q of order 
q, where c and d, in a less strict sense, correspond to the 
column and row weight of parity-check matrix, respectively. 
Based on this expression, we show that, when averaged on the 
whole ensemble, the fraction of codewords of small weight I 
in an LDPC code is at most asymptotically n _ r( c ~ 2 ) ( / 2 l as 
the coding length n goes to infinity. Next, using the upper- 
bound technique mentioned above, we analyze the asymptotic 
growth rate u gtC ,d(x) of the average weight distribution, where 
x denotes the normalized weight. A series of fundamental 
qualitative properties of u} q , c ,d[x) are found and proved. In 



particular, we show that for d > c > 3, 



J q.c 



d(x) has a 



unique zero xq in (0, 1 — l/q\. This zero just corresponds to 
the normalized minimum distance of a typical LDPC code, 
and hence provides important information about the code 
ensemble. Finally, we prove that for d > c > 3, there are at 
most a fraction Q(n~ r( c - 2 )V 2 l ) of all codes in the ensemble 
whose minimum distance is between the constant Iq and an, 
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where a £ (0, xq). 

The rest of this paper is organized as follows. In Sec- 
tion UH we introduce the notations and conventions to be 
used throughout the paper. In Section [TTU we define the 
ensemble of regular LDPC codes over a finite field and give 
its average weight distribution function; moreover, we study 
the asymptotic behavior of the average weight distribution for 
the small-weight case. The main analysis, consisting of two 
stages, for the asymptotic growth rate of the average weight 
distribution is performed in Sections HVl and IVl The minimum 
distance of individual codes in a regular LDPC code ensemble 
is analyzed in Section [VI] Section [VTTl concludes the paper. 

II. Notations and Conventions 

In this section, we introduce some basic notations and 
conventions to be used throughout the rest of this paper. 

• In general, symbols, real variables, and deterministic 
mappings are denoted by lowercase letters. Sets and 
random elements are denoted by capital letters. 

• The symbols Z, N, No, M. denote the ring of integers, the 
set of positive integers, the set of nonnegative integers, 
and the field of real numbers, respectively. For a prime 
power q > 2 the finite field of order q is denoted by ¥ q . 
The multiplicative subgroup of nonzero elements of ¥ q 
is denoted by F* . 

• The n-fold cartesian product of a set A is denoted by A n . 
An element of A n is denoted by x = (xi, X2, ■ ■ ■ , x n ), 
where Xi £ A denotes the ith component of x. 

• For any vector c £ F™, the weight w(c) of c is the number 
of nonzero symbols in it, that is, w(c) = \{i : ci ^ 0}|. 

• Given the functions / : X Y and g : Y — s- Z, their 
composite is the function g o f : X — > Z given by x H> 

<?(/(*))■ 

• Given the functions / : X\ — > Y\ and g : X 2 — > Y%, 
their cartesian product is the function fOg : Xi x Xi — > 
Yi x Y 2 given by (x 1 ,x 2 ) i-> (f(xi),g(x2))- 

• When performing probabilistic analysis, all objects of 
study are relative to a basic probability space (tt, A, P) 
where A is a cr-algebra in £1 and P is a probability 
measure on (il,A). For any event A g A, PA — P(A) 
is called the probability of A. Any measurable mapping 
of il into some measurable space (B. B) is generally 
called a random element. For any random set or function, 
we tacitly assume that their rt-fold cartesian products 
(e.g., A n or Oi=i F) are cartesian products of their 
independent copies. 

• All logarithms are taken to the natural base c and denoted 
by In. 

• For any x g [0, 1] and any integer q > 2, the entropy 
function H q (x) is defined by 

HJx) =xln- + (1- a;) In— \-xln(q - 1). 

x 1 — x 

For any x,y £ [0, 1], the information divergence function 
D(x\\y) is defined by 

D{x\\y) = x In — + (1 — x) In . 

y i - y 



• For any real functions f(n) and g(n) with n £ N, the 
asymptotic 0-notation f(n) = Q(g(n)) means that there 
exist positive constants c\ and c 2 such that 

cis(n) < f(n) < c 2 g(n). 

for sufficiently large n. 

• For x £ E, L^J denotes the largest integer not exceeding 
x, and \x~\ denotes the smallest integer not less than x. 

III. Regular LDPC Codes over Finite Fields 

We first define some basic F g -linear transformations. 

Definition 3.1: A single symbol repetition with a parameter 
c £ N is a mapping /^f p : F g — s- F^ given by x n- 
(x, x, . . . , x). 

Definition 3.2: A single symbol check with a parameter d £ 
N is a mapping : ¥ d q ->• F g given by x ^ £\ =1 Xl . 

Definition 3.3: A single symbol random multiplier map is a 
random mapping F q 7iM : ¥ q — >• ¥ q given by x — > Cx where C 
is an independent random variable uniformly distributed over 

Definition 3.4: A uniform random interleaver of F™ is a 
random automorphism E g n : F^ — > ¥ q given by x 4 
(xjj-im, Xn-xf2)i ■ ■ ■ > x Ji- x {n))i where II is an independent 
random permutation uniformly distributed over the symmetric 
group S„, i.e., all permutations on n letters. 

Next, we define a random linear transformation based on 
the above simple maps. 

Definition 3.5: F^ d n : F™ — > ¥ q n ^ d is a random mapping 
defined by 

p LD A fC HK pRM nV f REP fU 

r q,c,d,n ~ J q,d,cn/d r q.cn u ^q,cn u J q,c,n [ L > 

where c, d £ N, d divides cn, and 

n n n 

fKEP A /7N rREP fCHK ^ /— >. „CHK pRM A „rm 

J q,c,n \Ls 9' c ' •* q,d,n •> q,d s q-n vly Q 

i—1 i—1 i—1 

Considering the kernel of F^ d n , we thus obtain an ensem- 
ble of regular LDPC codes over ¥ q , which is called a random 
(c, d)-regular LDPC code over ¥ q and is denoted by C^? d Q 
This ensemble was originally introduced in O, fl2l . Ifl3ll by 
the method of bipartite graphs. 

To see the connection of F q ® d n with a bipartite graph, we 
may regard each f^f p as a variable node with c sockets and 
each j as a check node with d sockets. Then in total there 
are nc variable sockets and nc check sockets. We say that the 
ith variable socket and the jth check socket are connected 
by an edge if j = II (i), where II is the random permutation 
defined in Definition 13.41 We also define the label of the edge 
connecting these two sockets to be the random variable C 
defined in Definition 13.31 Then we dispose of the sockets (i.e. 
edges are considered as connections between variable nodes 
and check nodes). The resulting random graph (which may 
have repeated edges) is exactly the random regular bipartite 
graph with independent and uniformly distributed random edge 
labels taken from F^ as in (SI. 

'We shall tacitly assume throughout the paper that the block length n 
always takes values such that d divides cn. 
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Now let us investigate the weight distribution of C q n ^ d . The 
next theorem gives its average weight distribution. 

Theorem 3.6 (cf. ^j, ifTTV ): For c,d G N, the average 
weight distribution of C q ^ d is given by 



E 



dn) 
q,c,d 



(0 



(flcoeffry^s),^) 



(2) 



where ^ n c d (0 denotes the number of codewords of weight 2 
in C^? d (0 < I < n), coef (p(x),x l ) denotes the coefficient 
of x in the polynomial p(x), and 



&*) = ^ {I 1 + (9 " 1)*]" + (?-!)(!- a)"}' 
Furthermore, we have 



1 



^ d (0 <^(-)+ C /3 cn (c/) 



where 



w 9 ,c,d(a;) = H q (x) + - [S qi a(x) - lng] 



5 q Ax) = inf 5„ d {x,x) 

tfg,d(:Ei i) = rf£>(x||i) + p q ,d{x) 

qx 



P q ,d{x) = In 



1 + (?-!) 1 



9-1 



(3) 

(4) 

(5) 
(6) 
(7) 
(8) 

(9) 



Proof: The average weight distribution (O is in fact a 
known result. Note that 



E 



and 



{eeker/^K /d :K6)=d} 



1(9-1) 



^{cecg d | w (c) = /} 

For a proof of 

{c G ker /g c ™ /d : «,(£) = d}| = coef (^(z), ^ 

the reader is referred to JU Appendix III], Jg), ifTTl . 

Now let us prove the inequality (@). By the upper-bound 
technique introduced in Section |U it follows from (ffjl that 



E 



< 



)(q-l)( c - 1 »x cl 



for any x > 0. Taking 



(<?-l)(l-£) 
where i G (0, 1), we obtain 

(9-l)'(?)<fi 7 %) 



< 



(> d (i-) £ 



(10) 



where 



1 + (? - 1) 1 



qx 



Taking logarithms of both sides of ( TTOb and using the lower- 
bound in Lemma IA.1I we further have 
1 



\nE 



A 



(n) 
q,c,d V 



(/) < jff,(a) + - [*,,d(a,x) - ]nq\+cf3 cn (cl) 



where a = l/n. The theorem is finally established by taking 
the infimum of the right side over all x G (0, 1). □ 
Remark 3.7: Loosely speaking, for any a G [0,1], if we 
take I = an, then it follows from [4, Theorem 1] that 



lim — In coef ( g 

n— >oo ji \ ' 



(cn/d) 
q,d 



(x),x c 



— mi In — — -, — 

d x>o x da 

(crn/d) 

y q ,d 

gcmc 



= inf — In 

x>o m 



(,;) 



for any m > 0. Comparing this identity with the proof of 
Theorem [52] and noting that the second term in the right hand 
side of (|4]i is asymptotically negligible, we immediately have 



lim -InE 

n—>oo ji 



&q,c,d 



(«)■ 



The function cu qiCi d(x) thus represents the asymptotic growth 
rate of the average weight distribution of C q n } d , and hence 
deserves further investigations. In the subsequent sections, we 
shall provide an in-depth analysis of tJ q>Ct d(x), 

Although in general the average weight distribution of C q n ^ d 
is very complex, it becomes simple for some special d. The 
next two theorems give its complete characterization for d = 
1,2. 

Theorem 3.8: 

1 1 = 

otherwise. 

Proof: For d = 1, we have F^ d n = F™o£ 9 , c „o/ 9 R c % 
which is injective. In other words, the defining parity-check 



E 



A 



',i(0 



matrix of Cg C i has rank n, so that C qc ! ± = {0}. 
Theorem 3.9: 



E 



("Xci/2) 



d is even 
otherwise 



InE 



4?i,2(0] <(l-^)^(-)+c/3 c „(c0 



□ 



(11) 



(12) 



Proq/: By © it follows that 



9 v i(x) = \l + (q-l)x* 



Then we have 



coe{( 9q f 2 \x),x^)=l fe-^Cd^ 







d is even 
otherwise. 



This together with <j2j gives dTTT >. which further yields ([l2l i by 
Lemma IA.1I □ 
As shown above, the average weight distribution of C q n ^ d 
is trivial for d = 1,2. In the sequel, we shall therefore 
concentrate on the general case of d > 3. 
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Another well-known fact to be noted is that when q = 2 and 



d is even, the weight distribution of d satisfies A y ^ } c d (l) 



,(«) 



^4™id( n — f° r < / < n. This property simply follows 
from the fact that for even d the all-one vector is a codeword 

(n) 

of C 2 c d . In particular we have the following: 
Remark 3.10: For even d>2, 



E 



= E 



2,c,d 

U2,c,d( x ) = w 2,c,d(l - x) 



(13) 
(14) 



We close this section with a theorem on the asymptotic 
behavior of the average weight distribution for the small- 
weight case. 

Theorem 3.11: For d > 3 and constant weight I > 1, 



E 



c = 1 and I = 1 

9 = 2 and cl is odd 

6 (n _ r(c-2);/2l j otherwise. 



Proof: The trick of the proof is to find a precise ap- 
proximation of coef(g q °2 \ x )i x cl ) in (f2]l and to prove it by 
induction. For convenience, we define 

A(n,m) = coef (g^lix), : 

After some algebraic manipulations, we have 



e ;W)« 



where 



B(i) = 



(g-l)* + (-!)<(<?-!) 



Then it is observed that 

min{m,c£} 



A(n + l,m) = ^2 

i=0 

= A(n, m) 



A(n, m — i)B{i) 



min{m,c/} 



E 



A(n, m — i)B(i). 



Hence we have 

A(n,0) 
A(n,l) 
A(n,2) 



4(1,0) = 1 
4(1,1) = 

A(n-1,2) + 



A(n - 1,0)B(2) 



A(n- 1,2) 

efnLtr 



d(d-i)(?-i) 



A(n,3) 



A(n-1,3) + 
+ (^(n-l,fl)B(3) 



A(n -1,1)B(2) 




d(d-l)(d-2)(g-l)(g-2) 



6 



9 = 2 
otherwise. 



A(n, m) 



(15) 



We shall show by induction on m that 

q = 2 and m is odd 

otherwise. 

for all constant m > 2. Here, we only prove the general case 
of q > 2. The case of q = 2 can be proved by a similar 
argument with the fact B(i) = [1 + (— 1) 4 ]/2. Suppose that 
(fT31 > holds for 2 < m < k with fc > 3, then for m = k + 1, 

A(n, fc + 1) = A(n - l,fe + 1) 



E 

i=2 



A(n-l,fc-z + l)5(i) 



A(n - 1, fc + 1) + 6 ((n - 1) L( fc -!)/2J ) 



This asymptotic behavior implies that there exits a positive 
integer uq such that for n > no, 



A(n, k + 1) = A(n , k 



6 (r 




;L(fc-i)/2j 



L(fc+1)/2J 



Thus O holds for all m > 2. 

Finally, it follows from Theorem 13.61 and ( fT~5T > that 



E 



f)A{cn/d, cl) 




r(c-2)//2l> 



c = 1 and / = 1 
g = 2 and cl is odd 
otherwise 



□ 



as desired. 

Remark 3.12: The first and second cases of Theorem 13.1 II 
have the following alternative proofs: If c = 1 then the random 



code C q n ^ d , as the kernel of the reduced mapping f qd n/d ° 
F^f o Yj q:1l , has the same weight distribution as the kernel of 
fqd^n/d- ^ n Particular, C ™ c } d has no words of weight 1. If c 
is odd then every column of the parity-check matrix of C| j d 
(i.e. the transformation matrix of F^ d „) has odd weight. This 

implies that the all-one vector is in the dual code of cl" d and 
hence that all codewords have even weight. 



IV. Properties of the Function 5 qtd (x) 

As an important step towards understanding the function 
u: qcd (x), we analyze in this section the function 5 q ^{x) 
defined by ©. The proofs of lemmas in this section are 
presented in Appendix iDl 

In the sequel, we shall frequently use the following substi- 
tution to facilitate the analysis: 



1 



qx 



- z, qx 

2 = 1 -. (16) 

q- i 9-1 

Note that this transform is bijective and strictly decreasing, so 
we have 

a= (g-P(i-*), t= HM (17) 

9 9 
and z.z£ [—1/(9 — 1), 1] as x, x e [0, 1]. 
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Our first goal is to study the zeros of the partial derivative 
of S qt d(x,x) with respect to x. 

Lemma 4.1: For the function S qt d(x, x) defined by ©, 



dS q>d {x,x) _ dD(x\\x 



Ox 



where 



dpq,d(x) 

dx dx 
= qd{C, q , d {z) - z) 

(l-z)[l + (q-l)z\ 

a z + zJ- 1 + (q~2)z d 
\ + {q-l)z d 



(18) 
(19) 

(20) 



Lemma 14.11 shows that the zeros of dS qt d(x,x)/dx are 
determined by the equation Cq,d(z) — z = 0. We therefore 
proceed to analyze the function (q,d(z). The next three lemmas 
give the properties of C q ,d{z). 

Lemma 4.2: For q > 2 and d > 3, the function ( q .d(z) is 
continuously differentiable on [— l/(q— 1), 1] and its derivative 
is positive on (— l/(q — 1), 1). 

Lemma 4.3: For q > 2 and d > 1, 



CqX Z ) 



Z d -\1- z)[l + (q- l)z] 
l + (q-l)z d 



(21) 



Cq,d 



1 q — 2 and d is odd(22a) 



l 

9-1 



otherwise 



C g ,d(o) = o 

CU1) = L 



Lemma 4.4: Let 



- - 1 



g = 2 and <i is odd 
otherwise. 



(22b) 
(23) 
(24) 

(25a) 
(25b) 

zi{z) 



The equation £ qi d(z)—z — has a unique solution z\ 
in [— l/(q — 1), 1] for each z € [zi, 1] and has no solution in 
[— l/(q — 1),1] for z < Z\. The solution Z\{z) is continuous 
on [zi, 1] and is continuously differentiable on (z 1; l); its 
derivative is positive on (zi,l). Moreover, z\(z) E I q d( z )> 
where 



(z,0) 

{0} 

(0,z) 

{1} 



Equipped with Lemmas 14. 1144.41 we are now in a position 
to analyze the function 5 q d{x). 

Theorem 4.5: Let q > 2, d > 3, and 



z = Zl 

z € (zi,0) and d is odd 
z e (zi,0) and d is even 
z = 
ze (0,1) 
z = 1. 



q = 2 and d is odd (26a) 
otherwise. (26b) 



S q ,d{x) 



For the function 6 qt d(x) defined by ©, we have 

( In q x = (27a) 

X = 1 (27b) 

-oo xe (1— 3,1), <7 = 2, 

and d is odd (27c) 

ln(2d) -dH 2 (^J x=l-\,q = 2, 

and d is odd (27d) 

^ 5 q .d{x,£i) xe(0,xi) (27e) 

where p q ,d{x) is defined by (O and x\ = xi(x) is the unique 
root in (0,1) of the equation 

d5 q ,d(x,x) 



dx 



= o 



(28) 



solved for x as a function of x. The function x\{x) is continu- 
ously differentiable on (0,Xi) and its derivative is positive on 
(Q,Xi). Moreover, l\m x ^ + x\(x) — 0, lim^^- x\(x) = 1, 
and x\(x) £ I q> d{x), where 

f (3,1-1) IG (0> 1_1) 



Iq,d(x) 



qi v ' ql 

{1-1} X = l-1 

(x,l) 



x 6 (1 — ~, xi) and d is odd 
(1 — 1, x) x € (1— 1, cci) and d is even. 

The function <5 g ,d(x) is continuous on [0,Xi] and is continu- 
ously differentiable on (0,xi), in which case, 



dS q ^d(x) _ x(l - Si) 

dx xi(l — x) 

Proof: At first, Lemmas 07TJ 02] and 03] show that 

dS qid {0,x) 



(29) 



and 



<9x 

9^g,d(l,f) 



dx 



>0 V* 6(0,1) 
< Vx e (0,1). 



Therefore we have 



and 



8q,d(°) = lim ^9^(0, i) = p q ,d(&) 



(5g jd (l) = lim 5, )d (l,i) = p q ,d{l). 

X — ¥l~ 



This concludes (I27at and ( I27bb - 

A similar argument also shows that for odd d 



d5 2 ,d(x,x) 
dx 



< Vie 



1- 1 1 ) . , : ((). I ) 



so that 

fo.dOr) = lim 6 2 ,d(x,x) 



/ \ , l + (l-2x) d 

-dH 2 (x)+ hm In— ^— —1- 

£->l- (1 — X)" 1 - 1 x > 



= -dH 2 (x) +ln lim 



2(1 - 2x) d " 1 



-dH 2 (x) + In 



a-n- (l-x)(l-x) d ( 1 - ;r )- 1 
2 

(1 - x) lim^^i- (1 - 
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which yields J27c| > and j21d\ . 

For x € (0, x±), Lemma [4~4| shows that there is a unique 
zi = zi[z) G (-l/(<7 - 1),1) such that C,q,d(h) = z = 
l — qx/(q—l). Let x\ — (q—l)(l — zi)/q, which is essentially 
a function of x. Then it follows from Lemma [47X1 and l4~2l that 



and 



dS q ,d{x,x) 

dx 

dS q , d (x,x) 
dx 



< 



> 



Vx G (0,xi) 
V5 G (i i, 1). 



Therefore, 5 g ,d(a;) = xi), which concludes (127el >. 

Furthermore, Lemma l4~4l shows that xi(ir) is continuously dif- 
ferentiable on (0, x%) and its derivative is positive on (0, Xi). It 
also shows that \im x ^ a + x\{x) = and lim . - x%(x) = 1, 
and that &i(x) G I q ^d{x). 

Based on the above analysis, it is clear that S q ,d(x) is 
continuously differentiable on (0, x\). Furthermore, equation 
(127el > combined with Lemma |B~T1 gives (|29| i. 

Finally, let us show that 5 q ^d{x) is continuous at the end- 
points of the interval. Note that 5 q> d{x) is the infimum of a col- 
lection of continuous functions, so it is upper semi-continuous. 
Then it suffices to show that lim x _ > . + 5 q ,d{x) > S q> d(0) and 
li^-s-a- 5 qA x ) ^ S qi d(xi). Recall that lim a _ i . + xi(x) = 
and Viva. x ^ w - x\(x) = 1, so we have 

lim 6 q d(x) > lim p qd (xi(x)) = Inq 
lim S q ,d(x) > lim_ p 9)(J (i 1 (ar)) = /O g ,<j(l) 

and 

lim <52. d (x) > lim 



dH 2 (x) + In 
1 



l + (l-2xi(x)) c 
1 - xi(x) 



ln(2d) -dH 2 [- 



for odd d. The proof is complete. □ 
In Fig. Q] we give an illustration of the graphs of S q> d(x) for 
(q,d) = (2,5), (q,d) = (2,6), (q,d) = (3,5), and (g,d) = 
(3,6). 

V. Properties of the Function uj q:C , d (x) 

In this section, we proceed to analyze the properties of the 
function ui q , c ,d(x) defined by (|5j. Since LDPC codes are trivial 
when c > d, we shall sometimes assume c < d to exclude 
trivial cases. The proofs of lemmas in this section are presented 
in Appendix [E] 

At first, we calculate the value of ui q . c .d(x) at some special 
points. 

Lemma 5.1: Let q > 2, c > 1, and d > 3. 



w g ,c,d(0) = 0. 

1 i) (' 3) ln,/ - 



(30) 
(31) 



c c 

0V c ,d(l) = ln(g - 1) + -jPg.d(l) - 3 In g. (32) 
d a 



If g = 2 and d is odd then 
1 

d 

U)q,c,d(x) 



and 



^,ci(l-i) =(l-c)H 2 (] 



Vie l 



-lnd (33) 
d 



(34) 



Lemma 15.11 is an easy consequence of Theorem 14.51 so its 
proof is left to the reader. Next, let us calculate the first-order 
derivative of u)q tC ,d{x). 

Lemma 5.2: For the function u) qjCt d{x) defined by (0 with 
q > 2, c > 1, and d > 3, if x belongs to the case ( I27eb then 



du q ^d(x) 

dx 



= In 



1 -x 



1 — X\ 

xi 



+ M5-1) 
(35) 



which can be further expressed as 

duj q ^ d {x) _ J ! + (£?- l)5i 



In 



l-z? 



-i c-1 ' 



td-l 



(36) 

where ii is defined by d28l and 2i = 1 — qx\/{q — 1). 

The next lemma gives the value of du q c d(x)/dx at some 
special points. 

Lemma 5.3: Let q > 2, d > 3, and xi be defined by 




oo c = 1 


(37a) 


ln(d - 1) c = 2 


(37b) 


— oo c > 3. 


(37c) 


= 0. 


(38) 



If q = 2 and d is even then 

du q , c ,d{x) 



lim 

x — 



d.T 




If 5 ^ 2 or d is odd then 



lim 



dWq,c,d{x) 

dx 



c = 1 (39a) 
c = 2 (39b) 
c > 3. (39c) 



(40) 



To have more insights into Lj q<Ci d(x), we proceed to analyze 
the second-order derivative of u: qtC ,d{x). Since 

d 2 uo q , c ,d{ x ) _ dw q , c ,d{ x )\ dzi 

dzx \ 



dx 2 

and we note that 



di\ 
dx 



dx 
q dxi 



dx 



(41) 



q — 1 dx 

is negative on (0, a;i), our task is now to calculate the 
derivative d(dui q . c .d(x) / 'dx) / 'dz-y. 

Lemma 5.4: For the function ui qtCi d(x) defined by (0 with 
q > 2, c > 1, and d > 3, if x belongs to the case ( I27el i then 

d ( doj qiCi d{x )^ 
dzi \ dx 



(l-^)[l + (g-l)^[l + ( g -l)^ 



d-li 



(42) 
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i r 




~i 1 1 r 



q = 2,d = 5 
q = 2,d = 6 
q = 3, d = 5 
q = 3, d = 6 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 

Fig. 1. The graphs of & q , d {x) for (q, d) = (2, 5), (q, d) = (2, 6), (q, d) = (3, 5), and (q, d) = (3, 6) 



0.8 



0.9 



where 



^ c ,^) = £^[( c -i)(d-i)-i]z d - 



i=0 



-(g_l)[( C _l)(d-l)-l]i 
2d-3 

+ ( g -i)]Tz\ 

When c = 1, equation d42l reduces to 

d / du 2 ,cA x ) \ 1 



fi;j; 



(1 -£,)[! + 



(43) 



(44) 



We go on to analyze the function £ g c ,d(£) f° r 9 > 2, c > 2, 
and d > max{c, 3}. 

Lemma 5.5: For d > 3, the function £2,2, d(^) is positive on 
(— 1, 1). For q > 3 and d > 3, the function £ q .2,d(z) has a 
positive zero £2 in (—l/(q — 1), 1), and (, q .2,d(z) is positive 
on (— l/(q — 1), 22) and negative on (z 2 , !)■ 

For d > c > 3 with d even, the function £,2,c,d(z) has one 
zero z 2 in (0, 1) and the other zero z' 2 in (—1, 0), and £,2,c,d{z) 
is positive on (z' 2 , z 2 ) and negative on (—1, z' 2 ) U (z 2 , 1). 

For g > 2 and d > c > 3 with q ^ 2 or d odd, the function 
£ 9 ,c,d(2) has a positive zero z 2 in (-1/ (g-1), 1), and £q ;C ,d(z) 
is positive on {—l/{q— l),z 2 ) and negative on (£2,1)- 

We are now ready to give the qualitative properties of 



J q,c 



d(x). 



Theorem 5.6: Let q > 2, c > 1, d > max{c, 3}, and xi 
be defined by d26l i. The function w gjC) d(a;) defined by (O is 
continuous on [0,xi] and is twice differentiable on (0,Xi). 

If c = 1, then u>q, c ,d{%) is concave on (0,xi), and it is 
strictly increasing on (0,1 — 1/q) and strictly decreasing on 
(l-l/q, Xl ). 



If c = 2, then ui q:C: d(x) is strictly increasing on (0, 1 — 1/q) 
and strictly decreasing on (1 — 1/q, xi). Moreover, if q = 2, 
it is concave on (0, x\); otherwise, it is convex on (0, x 2 ) and 
concave on (x 2 , 1), where x 2 € (0, 1 — 1/q). 

If c > 3, q = 2, and d is even, then aj 9jC ( j(a;) is symmetric 
about the axis x = i. It is convex on (0,2:2) and concave 
on (#2, \ ) for some X2 S (0, i); it is strictly decreasing on 
(0,13) and strictly increasing on (x 3 , |), where xz € (0,^2); 
consequently, it has a unique zero x$ in (0, |], where xo G 
(X3, i], and it is negative on (0,Xo) and positive on (xq, 

For other cases, the function oj q , c ^{x) is convex on (0,^2) 
and concave on (x 2 , X\), where x 2 G (0, 1 — 1/q); it is strictly 
decreasing on (0,x 3 ) U (1 — l/q,xi) and strictly increasing 
on (#3, |), where x 3 g (0, X2); consequently, it has a unique 
zero xq in (0, 1 — 1/q], where xq € (^3, 1 — l/<?], and it is 
negative on (0, xq) and positive on (xq, 1 — 1/q). 

To provide an intuitive illustration of w q c d(x) in each case, 
the graphs of ui qtC d(x) for typical values of (q, c, d) are plotted 
in Figs. EHU 

Sketch of Proof: The proof is direct, and it depends 
on Remark 13.101 Theorem 14.51 Lemmas 15.1145.51 and identity 
(Put . Here, we only give the proof of the last paragraph of 
statements. 

Lemma 1531 and identity d4TT > show that uj q . c .d is convex on 
(0, X2) and concave on (x 2 ,x\), where x 2 £ (0, 1 — 1/q)- 
Furthermore, Lemmas 15.11 and 15.31 show that 







K) 









and 



,. duq,c,d(x) du qc ^ d {x) 
lim — —. = —00, 



>o+ dx 



dx 



= 0. 
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_q -j I I I I I I I I I I I 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Fig. 3. The graphs of u)2,cfi{x) for c = 1, c = 2, and c = 3. 
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Therefore, the derivative duj qyCy d(x)/dx has a unique zero 
X3 in (0,1 — 1/q), where £3 £ (0,2:2); it is negative on 
(0,13) U (1 — 1/q, x\) and positive on (x 3 ,l — 1/q). In 
other words, the function uj q ^ c ^{x) is strictly decreasing on 
(0, X3) U (1 — 1/q, x%) and strictly increasing on (x^, 1 — 1/q). 
The last statement about the unique zero in (0, 1 — 1/q) clearly 
follows. □ 
Remark 5.7: The zero xq in Theorem l5.6l iust corresponds 
to the normalized minimum distance of LDPC codes, in an 
average and asymptotic sense. It is in fact a function of q, c, 
and d, so we denote it by xo(q, c, d). We note that 

lim p q . d (x) =0 Vx € (0, 1) 
and hence for any r 6 (0, 1], 

lim x (q, \rd~] ,d) = x a , q . r 

where xo.g.r is the solution of H q (x) — r \nq = in (0, 1 — 
1/q). The detailed proof is left to the reader. Note that xo. q . r 
as well as the equation H q (x) — r In q = is closely related to 
the so-called asymptotic Gilbert- Varshamov (GV) bound over 
finite fields JT4J pp. 94-95]. This implies that regular LDPC 
codes with large c and d achieve the GV bound. 

VI. Minimum Distance of LDPC Codes 



Though we have shown in Remark 15.71 that regular LDPC 
code ensembles are asymptotically good, we are more inter- 
ested in the performance of individual codes of finite length. In 
this section, we shall investigate the minimum distance of an 
individual code in a regular LDPC code ensemble. To achieve 
this goal, we first establish an important inequality. 

Theorem 6.1: For q > 2, c > 1, d > 2, and x e (0, 1/q 2 ), 



W q , c ,d(z) < (J - l) XhlX + Kq,c,dX 



where 



hi(q- 



l) + -\n(d- 



3c. 



Proof: Put 



1 



(45) 



(46) 



(47) 



Then for any x £ (0, 1/q 2 ), x e (0, 1/q) C (0, 1 - 1/q). 
According to the definition (0 of Lj q>c ^d(x), we have 

u q>c ,d(x) 



< H q (x) + -{5 q4 {x,x) 
= H q {x) + cD(x\\x) 



Inq) 



In 



1 - 



1 - 



(a) 

< H q {x) + c£>(a;||x) + c 



qx 

q~i 

q(d 



2(9-1) 



-(c- l)H 2 (x)+x\n(q-l) 

, 1 / 1 m 1 , q(d-l) 

xm — + (1 — x) In — x + —. — : 

x 1 — x 2{q — 1) 

< (c — l)a; In x + x ln(q — 1) 



x A q(d-l) 
1-x 2{q-lY 
(c — l)x In x + x ln(g — 1) 



a; In — 

x 



a: In — 

x 



1 



2(9-1: 



(b) 



(c — l)ir In x + x ln(q — 1) 

1 , d- 1 x 
—x In 

2 x 



qx 



(d-l)(l-x) 2(q-l)_ 
< Q - l) x In x + (ln(q - 1) + | ln(d - 1) + 3c) x 



(c) 



where (a) follows from Lemma lAT2l and In a; < x—1, (b) from 
d47b . and (c) follows from q > 2, d > 2, and x < 1/q. □ 

Now, let us present the main result on the minimum distance 
of individual codes in a regular LDPC code ensemble. 

Theorem 6.2: For any code C C F™, we denote its mini- 
mum distance by d m i n (C). Then for q > 2, d > c > 3, lo > 1, 
and a G (0, 1 - 1/q), 

P{l < dmin(C q %) < na) 

< e L-r(c-2)0o+A)/2i j + q r% enUq , c , d (^ (48) 



where 



A = 



1 q = 2 and do is odd 
otherwise. 



(49) 



Proof: Since the minimum distance of a linear code is 
the minimum weight of its nonzero codewords, we have 



P{l < d min (C q %) < na] 

!\na\ 
(J [A ( q %{l) > l} 
l=lo 

\_7ia\ 

<E p {4"Uo>i} 



(b) 
< 



;=; 

io+3 



-r(c-2)J/21 



E 



e 



(nh n 



(c) 

< e 



(d) 

< e n 



-(c-2)(i +4)/2 > 



+ 6 n"e" u «^ w 



§ 9 r n -r(c-2)i D /2lj + (ni e nu q:C , d (^ 

where (a) follows from Markov's inequality, (b) from Theo- 
rems [32] and [3J~T] Lemma |A. II and the inequality l(n — I) < 
n 2 /4, (c) from Theorem 15.61 which shows that uj q . c .d(x) with 
x e [(Iq + A)/n,a] is upper bounded by either u; g . c .d((^o + 
4)/n) or Lu qyCy d(a), (d) from Theorem 16.11 and (e) follows 
from c > 3. 
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The above inequality holds in all cases. When q = 2 and 
d is odd, Theorem ECD shows that E[A^ d {l )} = 0, so 
we can further improve this inequality by simply replacing Iq 
with lo + 1. The proof is complete. □ 

Remark 6.3: If taking Iq = 1 in Theorem 16.21 we have 



P 



{<WO < no] 

< e (n-r( c - 2 )/ 2 i) + e (nh nu «>*^ E (50) 

Recall that uj q:C: d{x) has a unique zero xo(q,c,d) in (0,1 — 
!/<?), so we have 



P 



{d min (C^ d )<na}<e(n-r(- 2 )/ 2 l 



(51) 



for any a € (0, a;o(<7> c, <i)). Moreover, when c > 5, it 
follows from the Borel-Cantelli lemma that for any e > 0, 
the probability of the event 



1 



dmia(c£j jd ) < x (g, c,d) - e for infinitely many n 



is zero, so that 



' | lim inf — i 

n— >oo n 



(Cl d )>x (q,c,d) 



(52) 



The formula d50j >, for q = 2, was first proved (in a slightly 
stronger form for a different ensemble) by Gallager in |[T) . As 
for the general case of q > 2, Bennatan and Burshtein first 
showed in [ 8 1 that there exists some 7 > such that 



P 



{d min (0 < n 7 } < 6 (n 1 -/ 2 ) 



which is clearly weaker than (1511 . In [9|, Como and Fagnani 
proved a result similar to (IBTl l. 

Compared with previous results, the advantage of Theo- 
rem 16.21 is that we can use it to obtain results much better 
than (f50b by removing bad codes from the original ensemble. 
This viewpoint is formulated in the following theorem, which 
is an easy consequence of Theorem 16.21 

Theorem 6.4: Let q > 2, d > c > 3, lo > 2, and a £ 
(0, 1 - l/q). Let $„ : {All subspaces of F™} -> {0, 1} be a 
test function of linear codes such that for every linear code 
C, $„(C) - 1 implies d min (C7) > /„. If £[$ n (C^ d )] > 
Q(<p(n)) for some map 4>(n) : N —¥ [0, 1], then 



^{d min (0<na|$„(0 = l} 



< e 



l -r(c-2)(i +A)/21 



e 



ri2 e 



=,d(a) 



0(n) 



(53) 



where A is defined by ( |49l >. 
The proof is left to the reader. 

Remark 6.5: A simple test function can be defined by 
checking whether the parity-check matrix of a linear code 
contains all-zero columns. Then <fr n (C) — 1 if and only if 
the parity-check matrix of C contains no all-zero columns. It 

'When q = 2 and c is odd, we have a tighter upper bound 6(n 2 - c ) + 
6(n5e°" 2 - c ' d '"'). But for simplicity, we ignore this special case. 



is clear that $ n (C) = 1 is equivalent to d m j n (C) > 2, so it 
follows from ( BTT i that 



E 



*n(0] = > 2} = 6(1). 

Consequently, we have 

^{d^ n (O^H $ «(0 =i } 



<e 



+ e( 



c,d(«) 



(54) 



VII. Conclusion 



We provided a thorough analysis of the average weight dis- 
tributions of regular LDPC code ensembles over finite fields. 
The primary results are Theorems 13.111 14.51 15.61 and 16.11 
which are important for any analysis of regular LDPC codes 
based on the weight distribution. Furthermore, we proved a 
general result (Theorem 16.2b on the minimum distance of 
individual codes in a regular LDPC code ensemble, which 
includes all previous results as special cases. 
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Appendix A 
Some Useful Inequalities 

Lemma A.1: For any n £ N, define the function 

Pn{l) = H 2 (£j - 1 In Q VJ = 0, 1, . . . , n. 

Then 



< ff n (0 < — In ^ +e(rr x ) V0</<n 

2n \ n J 

and/3„(0)=/3„(n)=0. 

Sketch of Proof: Using Stirling's approximation: 

-j e A " Vn > 1 

where l/(12ra+ 1) < A„ < l/(12n). □ 
Lemma A.2: For all x £ [0, 1] and d £N, 

(l-x) d <l-dx+ d{d ~ l) x 2 . 



Proof: The inequality holds trivially for d = 1. Now 
suppose d > 2, then by Taylor's theorem, it follows that 

(1-^ = 1-^+^-^-^ 
for some y £ [0, x}. This thus concludes the proposition. □ 
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Appendix B 
Derivatives of H q (x), D(x\\x), and p q ,d{x) 

Lemma B.l: 

1 



dH q {x) 

dx 
dD{x\\x) 

dx 

dD{x\\x) 

dx 

dp q ,d{x) 

dx 



= In- 



= In 



x 

x(\ — x) 
x(l — x) 

x — X 



ln(g-l) 



x(l — x) 

qd(l- 



qx 



1 + (?-!)! 



qx 



(i ■ 



where p q ^d{x) is defined by ([EJ. 
The proof is left to the reader. 



Appendix C 
Descartes' Rule of Signs 

Theorem C.l (Descartes' Rule of Signs): If the terms of a 
univariate polynomial with real coefficients are ordered by 
ascending or descending variable exponent, then the number 
of positive roots of the polynomial (counted with their multi- 
plicities) is either equal to the number of sign changes between 
consecutive nonzero coefficients, or less than it by a multiple 
of 2. Since the negative roots of the polynomial equation 
f(x) — are positive roots of the equation f(—x) = 0, the 
rule can be readily applied to help count the negative roots as 
well. 

For a proof we refer the reader to lfT31 . 

Appendix D 
Proofs of Lemmas in SectionHVI 

Proof of Lemma \4.1\ By definition (0, ( TT~8T > follows 
immediately. Using Lemma IB. II and the change of variables 
( fl6l l yields 

qdz d - x 



dS q , d (x,x) 



qd(z — z) 



dx (1 - z)[l + (q - l)z] \ + {q-l)z d 

_ qd{z -z- z d - x - [(q -2)-{q - l)z\z d } 
{l-z)[l + {q-\)z][l + {q-l)z d ] 

qd{C q , d (z) - z) 



(l-z)[l + (q-l)z] 

as desired. □ 
Proof of Lemma \4.2\ To prove the lemma, we have 
to show that the derivative of Cq,d(z) is continuous on 
[— l/(q— 1), 1] and positive on (— l/{q— 1), 1). Some tedious 
manipulation yields 



CM 



[1 + (q - l)z d f 
where 

f{z) = l + (d- l)z d - 2 + (q- 2)dz d - 1 -(q- l)(d - l)z d 
~(q-l)z 2d - 2 . 



The continuity is obvious, even if q = 2 and d is odd. Our 
task is now to show that f(z) is positive on (— l/(q — 1), 1). 
The proof consists of two parts. 

First, we show that f(z) is positive on [0, 1). Note that the 
coefficients of f(z) have signs +, +,+,—, — . By Theorem lC.il 
it follows that f(z) has a unique positive zero. Since f(0) = 
1 > and /(l) = 0, it is clear that f(z) > for all I e [0, 1). 

Second, we show that f(z) is also positive on (— l/(q — 
1), 0) for both odd and even d. 

For odd d we have 

f(-z) = l-(d- l)z d - 2 + (q- 2)dz d ~ 1 

+ (q-l)(d-l)z d -(q-l)z 2d - 2 . 

If q > 3 then for all z, E (0, l/(q - 1)), 

f(-z) > \-{d-l)z d - 2 
(d-1) 



> 1 



> 



(q- l) d - 2 

2 d ~ 1 - {d - 1) 

(q- l) d ~ 2 



> 0. 

As for the case of q = 2, /(— z) reduces to 

1 - (d - l)z d - 2 + (d- l)z d - z 2d - 2 
which can be factorized as 

"d-3 



{1-zf 



E 

,i=0 



(i + !)(« + 2) 



'z* + z 2d - 5 - 1 ) 



(55) 



so that f(-z) > for all z £ (0, 1). 
For even d we have 



f(-z) = l + (d-l)z 



d-1 



(q-2)dz 



d-1 



(q-l){d-l)z d -(q-l)z 



zld-2 



> z d - 2 + (d 
d-1 



If 



~d-2 (q-2)d^ d - 



-z 



d-1 



_ s d-2 



= 



□ 



for all z € (0, l/(q — 1)). The proof is complete. 

Sketch of Proof of Lemma \4.3\ Identity (ETT i is proved 
by a straightforward argument using definition (1201 1. Equations 
( I22bb . (1231 1. and (l24l i are immediate consequence of (BTT i. As 
for (I22al i. we note that (1211 with q = 2 and odd d gives 



v d-l 



C2,d(-1) + 1 = 



(1 



1-Z 



+ z 



d-1 



2 
d 



so that C2,d(-1) = 2/d - 1. □ 
Proof of Lemma \4.4\ Lemmas 14.21 and 14.31 show that the 
range of C«,d(2) for z £ [-!/(? - 1), 1] is 



g-i 



= [^1,1] 



and therefore the equation C2.d(^) — z = has a unique 
solution in \—l/(q — 1),1] for each z G [zi,l] and has no 
solution in [— l/(q — 1), 1] for z < z\. 
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Since ( q ,d(z) is continuously differentiable on [— l/(g — 
1), 1] and its derivative is positive on (— l/(g— 1), 1), it follows 
from the inverse function theorem that the solution z\{z) is 
continuously differentiable on (z%, 1) and its derivative is also 
positive on (z\, 1). The continuity of Zi(z) at endpoints also 
follows. Moreover, Lemma 1431 shows that 



By the symmetric property (Remark |3.10l ), we also obtain (f39b - 
From Theorem 14.51 it follows that 



q) g-i 



if z 



(q,d( Z ) = Z l 
CqA Z ) > Z 
(q,d( Z ) < Z 

C*,d(*) = o 

CqX Z ) > Z 
CqA Z ) = 1 

This implies that z\{z) S I' q d (z) 



l 

9-1 



if z G (— ^ry,0) and d is odd 
if z € (— ^zj,0) and d is even 
if z = 
if ze (0,1) 
if z = 1. 



This together with equation ( l36l l gives 

Again by Theorem 14.51 it follows that lim - zi = 
-l/(g - 1). Then d36j with g ^ 2 or d odd gives gQj. □ 
Proof of Lemma 15.41 It follows from Lemma 15.21 that 

d ^ duj q ^ d (x) 



dx 



q 



q{c-l){d-l)zt 



d-2 



□ 



Appendix E 
Proofs of Lemmas in SectionM 

Proof of Lemma \5.2\ Definition Q and equation ( f29b 
show that 



duJq, c ,d( X ) _ dH q (x) 



dx 



= In 



dx 

1-x 



■In 



x(l — ii) 
Xi(l — x) 



ln(g - 1) 

e-1 



• 111 



x(l — X\) 
Xi(l — x) 



1 — X\ 
Xi 



[l + (q- l)zi](l - h) (1 - ^-^[l + (q- I)**" 1 ] 

= g^dCg) 

(l-^Il + fe-lJzJIl + Cg-l)^- 1 ]' 
This concludes (142b . while the first equality with c = 1 gives 

61. □ 

Proof of Lemma [Or Since £g lC) d(0) = 1, it suffices to 
determine all zeros of £g, c ,rf(z) in (— l/(g — 1), 1). The proof 
consists of two parts. 

First, we check the zeros of £q tC ,d( z ) in (0, 1). We note 
that the coefficients of £, q , c ,d{z) have signs +,...,+,—, 
— ,+,...,+. By Theorem IC.ll it follows that £ q , c .d( z ) has 
zero or two positive zeros. On the other hand, 



where (a) follows from Lemma |B~T1 By Lemma |4~T1 equation 
( |28T i is equivalent to ( q ,d(zi) — z — 0, where z — 1 — gx/ (q — 1) 
and li = 1 — q£i/(q — 1). After some manipulations, we obtain 



and 



&a,d(l)= o(9-2)(rf-l)(rf-2). 



x 
xi 



1 - z 



d-1 



i-2i l + (g - 



and 



1 



l + l + (g-l)i 



rf-i 



1 — Xi 



Then 

<fc; g ,c,d(a:) 
dx 



In 



In 



l + (g-l)£i 

f (g-l)(l-ai: 
I ii 

1 + (g - 1)Z! 



1 + (g - 



x(l — Xl) 



1-Z! 



Xx(l 
1 - 



zt x 



l + (g-l)z 



d-l 



The proof is complete. □ 
Proof of Lemma 15.51 From Theorem 14.51 it follows that 
lim a: _ ) .o+ z i — 1- Then equation d36l > with c = 1 and c > 3 
gives d37ab and (I37ct . respectively. As for c = 2, we have 

dw g , Cj d(a;) 



lim 

x^o+ dx 

= lim In 

21— >i~ 



lim In 

Zl— vi— 



[l + (g _ 1)^(1-^) 

(l-^)[i + (g- l)sf x \ 



l + (g-l)i 



d-l 




ln(d - 1) 



Then for q > 2 and d > c > 3, £ g , c ,d(<2) has a unique zero 
£2 in (0, 1). As for c = 2, £5,2,^(2) with g > 3 has a unique 
zero z 2 in (0, 1) since £g,2,d(l) = and ^ 2 d (l) > 0, while 
&,2,d(z) has only one zero z, = 1 in (0, oo) since £ 2 2 d (l) = 
(a zero of multiplicity 2), so that %2,2,d(z) is positive on (0, 1). 

Second, we check the zeros of £ q . c ,d( z ) in (— l/(g— 1), 1). 
To facilitate the analysis, we consider the function 

fg,c,d( z ) = (! + «)£g,c,d(-2)- 

Then the zeros of £, q . c .d(z) in (— l/(g— 1), 0) are just the zeros 
of f q . c . d (z) in (0,1/(5-1)). 
If d is odd, we have 

/,,c^) = (l-^- 1 )[l + (?-l)^- 1 ] 

+ (c - l)(d - l)z d - 2 (l + z)[l - (g - l)z] 

which is clearly positive for all z E (0, l/(g— 1)). 
If d is even, we have 

/ g , c ,^) = (l + i d - 1 )[l-(5-l)2 d - 1 ] 

- (c- l)(d- l)z d ' 2 (l + i)[l - (g - 1)2} 
= 1 - (c- l)(d- \)z d - 2 

+ (g-2)[( C -l)(rf-l)-l]z d - 1 

+ (g - l)(c - l)(d - l)z d - (g - l)z 2d - 2 . 

When g = 2, it reduces to 

h,cA z ) = l-{c-l)(d-l)z d - 2 +(c~l){d-l)z d -z 2d - 2 . 
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Since the coefficients of f2, c ,d(z) have signs +,—,+,—, it 
follows from Theorem IC.ll that f2,c,d{z) has one or three 
positive zeros. Moreover, we note that 



and 



f 2 ,c,d(°) = !. hcA 1 ) = °> /2,c,«l(oo) = 



£, c , d (l) = 2(c-2)(d-l). 



Then f2, c ,d{z) with c > 3 has a unique zero z' 2 in (0, 1), while 
f2,2,d{z) is positive on (0, 1) because of d55b . Finally, let us 
show that fq tC ,d(z) is positive on (0, l/(q — 1)) for q > 3, 
c > 2, and d > max{c, 4}. Since q > 3, c > 2, d > 4, and 
5 < 1/(9-1), 



f q ,cd{z) > 1 - (c - l)(d - l)z d - 2 (l - £ - 2z 2 ) 



- z d - 1 - z 2d - 3 



> 1 - (c- l)(d - 1)5 



d-2 



(56) 
(57) 



For z G (0, |], inequality d57| i shows that 

UcAz) > l-(d-l) 



For z€ (i |], inequality (TS6b shows that 



fq,c,d(z) > 1 



4(d- l) 2 /2 



>l--. 3 ' 
893 



d-2 



2\ d— i ^ 



3125 



For z € (| , i), inequality d56l l shows that 



f q , c ,d{z) > 1 



7(d- l) 2 (l\ d ~ 2 flV' 1 /1 N 2 '' ; 



25 



7 



^ 1 — • 3 



25 V 2 

- Ill 
~ 800' 

The proof is complete 



2 j \ 2 

■'< /i\ 5 

2 



□ 
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